Aging and Immortality in a Cell Proliferation Model 
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We investigate a model of cell division in which the length of telomeres within a cell regulates 
its proliferative potential. At each division, telomeres undergo a systematic length decrease as 
well as a superimposed fluctuation due to exchange of telomere DNA between the two daughter 
cells. A cell becomes senescent when one or more of its telomeres become shorter than a critical 
length. We map this telomere dynamics onto a biased branching diffusion process with an absorbing 
boundary condition whenever any telomere reaches the critical length. Using first-passage ideas, 
we find a phase transition between finite lifetime and immortality (infinite proliferation) of the cell 
population as a function of the infiuence of telomere shortening, fiuctuations, and cell division. 

PACS numbers: 02.50.-r, 05.40.Jc, 87.17.Ee 



I. INTRODUCTION 

Aging is a complex and incompletely understood 
process characterized by deteriorating cellular, organ, 
and system function. Replicative senescence, the phe- 
nomenon whereby normal somatic cells show a finite pro- 
liferative capacity, is thought to be a major contributor to 
this decline [ij. Telomeres are repetitive DNA sequences 
((TTAGGG)„ in human cells) at both ends of each lin- 
ear chromosome and their role is to protect the coding 
part of the DNA. Normal human somatic cells become 
senescent after a finite number of doublings. 

As a general rule, cells for which telomerase activity 
is absent lose of the order of 100 base pairs of telomeric 
DNA from chromosome ends in every cell division. This 
basal loss has been attributed to the end replication prob- 
lem '5] in which the DNA-polymerase cannot replicate 
all the way to the end of the chromosome during DNA 
lagging strand synthesis. Additional post-replicative pro- 
cessing of the telomeric DNA is necessary to protect the 
end of the chromosome from being recognized as a dou- 
ble strand break in need of repair. This processing also 
contributes to the end replication problem. 

What remains unexplained, however, is why senescent 
cells occur in cell cultures long before the expected num- 
ber of cell divisions estimated from the gradual basal loss. 
It has been shown recently that in addition to the basal 
loss of ~ 100 base pairs per division, a complex set of 
events that leads to telomere exchange between sister 
chromatids can occur. This telomere sister chromatid ex- 
change (T-SCE), together with basal telomere loss and 
a number of observed or suggested telomere recombina- 
tion events, collectively define telomere dynamics, and 
this dynamics leads to a wide distribution of telomere 
lengths in cell cultures [3]. Telomere exchange can also 
occur between the telomeres of different chromosomes. 
Currently available data cannot distinguish between this 
process and T-SCE. It is also believed that sister chro- 
mosome exchange is induced by DNA damage. Because 
the sister chromatids are in closer proximity compared 



to the distance between different chromatids, it can be 
hypothesize that the probability for telomere interchro- 
matid exchange (T-ICE) is smaller than the probability 
for sister chromatid exchange between sister chromatids. 

Recently, one of the present authors proposed a theory 
[1], based on telomere dynamics including T-SCE and 
T-ICE, that is capable of explaining Werner's syndrome, 
an inherited disease characterized by premature aging 
and death, and a subset of cancers that seem to use a 
recombination mechanism to maintain telomere length. 
This latter mechanism, known as alternative lengthen- 
ing of telomeres (ALT), as well as additional telomerase 
activity in some cases, are thought to contribute to the 
large proliferative potential of these cells ^5] . Here prolif- 
erative potential denotes the time over which cells can 
continue to divide. In Ref. we showed that both 
Werner's syndrome and ALT can be described within the 
general framework of telomere dynamics in which there 
is an elevated rate of exchange of telomere DNA between 
two daughter cells, as observed in many experiments 
One of the main numerical results from Ref. Q is a tran- 
sition from finite to infinite proliferative potential in a 
cell culture as the parameters controlling the telomere 
recombination rates are varied. 
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FIG. 1: Schematic illustration of telomere evolution. The 
telomeres in the initial cell contain a; = 15 units. Upon di- 
vision, each telomere ostensibly shortens by Ax — 5 units, 
but additional exchange of 5x = ±2 units between daughter 
telomeres (dashed) leads to final lengths of 8 and 12. 
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In this work, we analytically investigate an idealized 
model of cell proliferation in which the number of cell 
divisions before senescence occurs is controlled by the 
dynamics of telomeres during cell division. Each cell con- 
tains a certain number of telomeres. When a cell divides, 
the telomeres in each daughter cell ostensibly shorten by 
a fixed amount Ax. In addition to this systematic telom- 
ere shortening, the effect of T-SCE processes during cell 
division leads to a superimposed stochastic component to 
the telomere length dynamics by an amount 6x (Fig. [1]). 
Thus in each cell division event, the length of an indi- 
vidual telomere evolves by the combined effects of these 
systematic and stochastic mechanisms. 

When the length x of any of the telomeres within a cell 
reaches zero, the cell stops dividing and becomes senes- 
cent. On the other hand, the stochastic component of the 
telomere dynamics provides the possibility for a telom- 
ere to occasionally grow when a cell divides. This sub- 
population of cells with long telomeres and thus higher 
proliferative potential can become even more so at the 
next cell division, a mechanism that allows a long-lived 
subpopulation to thrive. 

We are interested in basic statistical properties of the 
cell proliferative potential. Some fundamental questions 
that we will study include: 

1. Can a cell population divide indefinitely? 

2. How long does it take for a cell population to be- 
come senescent? 

3. How many dividing cells exist after a given number 
of divisions? 

Within an idealized model of telomere kinetics described 
by Eq. ^ for a single telomere per cell, it is possible 
to answer these questions analytically by mapping the 
telomere dynamics to a first-passage process. Using this 
approach, we find a phase transition between a finite- 
lifetime cell population and immortality as a function 
of three basic control parameters — the magnitude of the 
systematic part of the telomere evolution, determined by 
Ax, the effective diffusion coefficient associated with the 
stochastic part of the telomere evolution, determined by 
{6x{t)'^), and the cell division rate. 



II. TELOMERE REPLICATION MODEL 

In our telomere replication model, we assume that the 
initial length of each telomere in a cell is Xq. In any cell 
division event, the length of a telomere changes by two 
distinct processes: 

i. a systematic shortening of each telomere by Ax; 

ii. an additional stochastic component of the length 
change of magnitude Sx. 

Thus the length of a telomere changes according to 

X — !■ X — Ax + Sx. (1) 



Here the stochastic variable Sx accounts for the T-SCE 
processes that we assume to have mean value equal to 
zero, (Sx) = 0, and no correlations at different times, 
{Sx{t)Sx{t')) cx S{t — t'). The justification for the absence 
of correlation is that T-SCE events have been linked to 
DNA damage [7|, which occurs randomly in the cell. 

Because of the systematic and stochastic contributions 
to the change in telomere length in each division event, 
the length of a telomere undergoes a biased random walk, 
with a bias toward shrinking. It is instructive to estimate 
the relative importance of the systematic and stochastic 
components of this length evolution. For this purpose, 
we define the time unit as the physical time between cell 
divisions St. Thus in the absence of stochasticity a telom- 
ere shrinks to zero length in xo/Ax cell division events. 

It is now helpful to recall some basic numbers about 
human telomeres to connect our mechanistic telomere 
model and real cell division: 



parameter 


definition 


numerical value 


St 


time between 
cell divisions 


20 minutes ~ 
several hours 


Xo 


initial telomere 
length 


» 10* base pairs 


Ax 


systematic length 
decrease per division 


« 10^ base pairs 


Sx 


stochastic length 
decrease per division 


« 10"^ base pairs 



The quantity Sx is known only very roughly. It is possi- 
ble, with low probability, that even whole telomeres can 
be lost in TCE processes Q. 

With the above numbers, a telomere shrinks to zero in 
xo/Ax « 10^ = N cell division events with purely deter- 
ministic shrinking. Now consider the role of stochasticity: 
in N cell divisions, the root-mean-square length change 
due to stochastic events is ^rms = \/ N{Sx)^ « 10'^. Thus 
in the time for a telomere to systematically shrink from 
lO** to zero, stochasticity gives a length uncertainty of 
10'^ — a 10% correction to the bias. 

Finally, we need to include the role of cell division on 
this biased random walk description to arrive at a theory 
for telomere dynamics. That is, we need to allow a ran- 
dom walk to replicate as it undergoes biased hopping. 
While there are many ways to parameterize the effects 
of bias, stochasticity, and replication in the continuum 
limit, all such models lead to the following convection- 
diffusion equation with multiplicative growth: 



dn(x, t) 
dt 



kn{x, t) + V 



dn{x,t) d^n{x,t) 



dx 



+ D 



dx"^ 



(2) 



in which v represents the bias for telomere shrinking, D 
accounts for the stochastic part of the telomere length 
evolution, and k accounts for cell division. While there 
is only an indirect connection between the model param- 
eters v^D^k, and the parameters Ax,(5x that account 
for what happens to a telomere in a single cell division, 
this continuum description has the advantage of captur- 
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ing the physical essence of telomere dynamics while being 
analytically tractable. 

The basic question that we seek to understand is how 
long it takes for a cell to become senescent, an event that 
occurs when the length of one of its telomeres reaches 
zero. This condition translates to an absorbing bound- 
ary condition at a; = for the biased diffusion process 
that described the telomere length distribution. We now 
exploit some classic results about the first-passage prob- 
ability of biased diffusion (see Appendix and [9]) to de- 
termine the evolution of the telomere length distribution. 
Let the initial number of cells be Nq. The solution of ([2|) 
is simply 



(3) 



where c{x,t), given in (|A.2p . is the solution to the 
convection-diffusion equation with an absorbing bound- 
ary condition at a: = and the initial condition that all 
cells have initial telomere length xq. 

We may also treat the situation in which each cell con- 
tains M independent telomeres as a zeroth-order descrip- 
tion for cells that contain many telomeres whose dynam- 
ics is coupled by the T-SCE exchange process (schemat- 
ically illustrated in Fig. [T|) . This assumption of indepen- 
dence of different telomeres allows us to apply the single 
telomere per cell first-passage description with only mi- 
nor modifications. For cells containing M independent 
telomeres of lengths x = {xi, . . . ,xm), the density of cells 
n(x, t) satisfies the M-dimensional convection-diffusion- 
growth equation 



9n(x, t) 
dt 



fcn(x, t) + wVn(x, t) + DV'^nijn, t) , (4) 



with absorbing boundaries when any telomere length Xi 
reaches 0. The solution simply factorizes as a product of 
one-dimensional solutions 



M 



n{-K,t)^NQe^*Wc{x„t) . 



(5) 



1=1 



If each telomere has initial length xq, all the c{xi^t) are 
the same and are given by (|A.2|) . We can also straightfor- 
wardly study the case where each telomere has a different 
initial length by merely using the unique initial length of 
each telomere in Eq. ([5]). It is worth mentioning that 
in addition to its application to cell division statistics, 
Eg. (|3]) is also related to the Fleming- Viot (FV) process 
in which a population of diffusing particles can get 
absorbed at boundary point and then be re-injected into 
the system at a rate that is proportional to local particle 
density. 

From our description of telomere dynamics as a biased 
branching-diffusion process, we now determine basic fea- 
tures about the time dependence of the cell population 
and the statistics of telomere lengths. 



A. Number of Dividing and Senescent Cells 

Cells in which each telomere has positive length can 
divide. The number of such active cells is given by the 
integral of the number density of cells over the positive 
2*^-tant {xi,X2t ■ ■ ,xn > 0) of length space: 

A^active(i) = / «(x, t) dx - Noc'"' S^' {t) , (6) 
Jx>0 

where S{t) is the survival probability of a biased ran- 
dom walk (given by (IA.4|) '). From (jA.6p . the long-time 
behavior of Afactivc(i) is given by 



N., 



active 



it) 



f .1 ^ ^^0 '^^^Mvxa/2D ^-3A//2 ^{k-Mv'^ /4D)t 
\ \ IT v'^ J ' 

(7) 



oc t 



-3M/2 ^kt{<iM-l)/<^M 



with eM = ADk/Mv^ . Thus the fundamental parame- 
ters of the system are the Peclet number Pe = vxq/2D 
a dimensionless measure of the relative importance 
of the bias and the stochasticity and, more importantly, 
the dimensionless growth rate, cm- Note that both the 
diffusion coefficient D and the bias velocity v are propor- 
tional to fc, since the telomere length changes occur only 
when a cell divides. As a result, the growth rate em is 
actually independent of k. For ea/ < 1, cell division is 
insufficient to overcome the effect of inexorable death due 
to the systematic component of the telomere shortening 
and the population of dividing cells decays exponentially 
in time. 




100 



FIG. 2: Plot of the number of active cells versus time from 
the second line of Eq. ((Tj for the case of M = 1 telomere per 
ceU, with ei = 0.8 (solid) 0.9 (dashed), and 0.99 (dotted). 



To get a feeling for the interplay between telomere 
shortening and cell division, let us employ the numeri- 
cal parameters given at the beginning of this section in 
Eq. ([7]) to obtain the total number of dividing cells. Since 
cells double at each time step, A: = In 2 when we express 
t in units of cell division times. Furthermore, we define 
the coefficients a and /3 by ti = a x 10^ and D = f3 x 10*. 
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Then the time dependent exponential factor in ([7]) for 
the case M = 1 becomes 



,(ln 2-0^4/3)* 



The exponent can be either positive or negative exponent 
depending on a and /3, which, in turn, depend on de- 
tails of the telomere evolution a single cell division event. 
Thus, using the numbers appropriate for humans, senes- 
cence or immortality is controlled by details of telomere 
evolution as encoded by the coefficients a and /3. 

We also obtain the total number of cells that become 
senescent at any given time during the evolution as the 
diffusive flux to Xi = for any i. For the symmetric ini- 
tial condition, the number of cells reaching any boundary 
Xi = 0, i — 1,2, . . . , M), are the same. Hence we may 
consider a single boundary, say i = \. The number of dy- 
ing cells at this boundary can be written as the integral 
over the a;i = surface 



Ji(i) 



>0 
I =0 



9n(x,t) 



oxi 



(8) 



where the integral is over all the M — 1 coordinates per- 
pendicular to xi. Note again the absence of a convective 
term in this expression because there is no convective 
flux when the concentration is zero. Using the product 
form ^ of r7,(x, t) for the total number of dying cells we 
obtain 



J{t) = MJi{t) = NoMe^'F{t)S'''-'{t) 



(9) 



where S{t) and F{t) are given by (|A.4|) and (|A.3p re- 
spectively. The total number of senescent cells that are 
produced during the course of the evolution is 



dt 



oo poo fi^^'^(i') 

J{t)dt^-No / dte'' ^ > 
Jo 

1-k dte^' S^Ut) 
Jo 



(10) 



where we again use the fact that F{t) — —dS{t)/dt to 
perform the integration by parts. 

For the special case of one telomere per cell [M = 1) 
we have 

J(t) = e'=* g-(.o-„t)V4Dt_ 

In this case the total number of senescent cells is 

NqXo 



)^ f _J_ kt (xo-vtf/iDt 1. (.^^ 



\/4 

We now make the substitution z = to recast the 

above integral into the form 



A^s, 



^0^0 ^vxo/2D 



dz e 



^xlz^/AD-{fj,-k)/z' 



(13) 




FIG. 3: Plot of the total number of senescent cells that are 
produced during the evolution of the cell population as a func- 
tion of the dimensionless growth rate ei — ADk/v^. 



Usingthe ^--az^-b/z'- = 
Ref. [TjI, we thus obtain 



,-\/4ah 



from 3.325 of 



A^s, 



Noe^" exp 



cxp 



(14) 



which is plotted in Fig.[3]for the case A'o = 1 and Pe — 10. 

This result for Ascn holds only for ei < I; this is 
the regime where the cell population eventually becomes 
senescent so that the total number of senescent cells pro- 
duced during the evolution is finite. When ei ^ 1, 
the leading behaviors of the two exponential factors in 
Eq. dH]) cancel and Agen ^ (l + ei/2Pe). Thus as ei 
(no division), the initial cell immediately becomes senes- 
cent so that the total number of senescent cells that are 
produced equals one. This unrealistic result arises be- 
cause in the continuum description telomeres can shrink 
to zero length before any cell division can occur. On the 
other hand, for ei > I, the population is immortal and 
an infinite number of senescent cells are produced during 
the evolution. 



B. Telomere Length Distribution 

A curious aspect of the population of dividing cells 
is the dependence of the cell density on telomere length 
as t —^ oo. For any number of telomeres per cell M, 
the telomere length distribution is independent of this 
number, since the distribution is simply proportional to 
c{x,t) in Eq. (|A.2p . In the t — > oo limit, we then obtain 



c{x, t — > oo) oc 



XXq 



oc xe 



vx/2D 



-(v^ /AD)t-v(x-xa)/2D 



(15) 
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Thus apart from an overall time-dependent factor, the 
density of dividing cells in which the constituent telom- 
ere has length x is linear in x for small lengths and has 
an exponential cutoff for large lengths. From the distri- 
bution given in Eq. (|15[) . the mean telomere length goes 
to a constant for large times 



X c{x, t) dx AD 
c{x, t) dx V 



(16) 



independent of whether the total cell population is grow- 
ing or decaying. The variance of the telomere length also 
approaches a constant 



{x') - {xf 



8D^ 



(17) 



A more pedantic way to arrive at this same result is to 
compute the exact time-dependent mean telomere length 
{x{t)) using Eq. (|A.2p for c{x, t) and then taking the t 
oo limit. In summary, the mean telomere length is time 
independent, as predicted by the asymptotic form of the 
telomere length distribution in Eq. (fT5|) . We emphasize 
that this result pertains to the fraction of cells that are 
active. If all cells — senescent and active — are included in 
the average, then (x) would asymptotically decay with 
time. 



where 7 does not necessarily have to be an integer. Using 
this general formula the variance of the lifetime is 



2Dx() 



t;3(l-ei)3/2 



(20) 



In fact, the n^^ cumulant can be obtained simply from 
the general formula 13] 



Therefore for the higher moments of the lifetime, the dif- 
fusion coefficient does play an essential role. 

For M = 1, we immediately obtain from Eq. ((T4|) 



1 



(21) 



Thus the mean cell proliferative potential diverges to in- 
finity as ei — !■ 1 from below. This result for the number 
of cell divisions before senescence is one of our primary 
results. Notice that in the case of no cell division {k — 0, 
ei = 0) the average lifetime (t) — xq/v. That is, (t) 
coincides with the time for a biased diffusing particle to 
be convected to the origin. It is surprising at first sight 
that diffusion plays no role in determining the average 
number of cell divisions. Exactly the same type of result 
arises for the discrete random walk with a bias v 11411. 



C. Mean Proliferative Potential and Immortality 



D. Senescence Plateau 



As discussed above, for eA/ < 1, all cells eventually 
become senescent. However, for eM > 1, a subpopulation 
of infinitely dividing cells arise so that the average cell 
population becomes immortal. We make this statement 
more precise by computing the average lifetime of the 
population for < 1- This lifetime is defined as the 
average age of each cell when it becomes senescent and 
is thus given by 



tJ{t)dt / te^^F{t)S^'-^{t)dt 
° - ^0 , (18) 



J{t) dt 



e^'F{t)S'"-'{t) dt 



where we use (|9]) for J{t), the number of cells that become 
senescent at time t. More simply, the mean proliferative 
potential (lifetime) can also be obtained from the total 
number of senescent cells via (t) = ^ In A^son- 

Using the general definition of Eq. (fT5|) , it is straight- 
forward to compute higher moments of the lifetime for 
the case M = 1 . Here the integrals in (fTS]) can be written 
in terms of the modified Bessel functions of the second 
kind (see #12 in 3.471 of Ref. (H) to give 



7/2 



(19) 



A useful characterization of the number of cell divisions 
distribution function of a population is the senescence 
rate m(t). The senescence (or mortality) rate is the ratio 
of the number of cells that become senescent at time t 
to the total number of cells that are still dividing at this 
time. Equivalently, the senescence rate is the probability 
that a randomly chosen dividing cell becomes senescent 
at the next moment. The senescence rate is thus given 

by 



m{t) 



N{t) 



(22) 



where the number of dying cells J{t) and dividing cells 
N{t) are given by Eqs. (O and ^ respectively. Substi- 
tuting these expressions into p2l) , we find that the senes- 
cence rate is independent of the number of telomeres per 
cell M, i.e. 



m{t) 



F{t) d\nS{t) 



S{t) 



dt 



(23) 



Using the asymptotic form of S{t) given in (jA.6p . the 
senescence rate approaches a time-independent value in 
the long-time limit and is given by 



(24) 
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Amazingly, as the cell population ages, the senescence 
rate of the cells that remain dividing ultimately tends 
to a constant value for large times. This phenomenon is 
known as the mortality plateau. Namely, the probability 
that the somatic cells of an organism become senescent 
becomes independent of its age in the long-time limit. 
This surprising fact was observed experimentally in hu- 
man populations [1^1 and for fruit flies 16]. It was also 
observed numerically in a model of aging that is similar 
to ours [13] ■ The existence of such a senescence plateau 
is actually typical of a wide range of Markov processes 

M. 



III. SUMMARY 

We studied an idealized model for the dynamics of 
telomere lengths during cell division that is based on a 
systematic basal loss and a stochastic component to the 
evolution that arises from telomere sister chromatid ex- 
change. This model captures essential features of telom- 
ere dynamics in cell cultures. Because of the compet- 
ing influences of cell division, which obviously increases 
the number of proliferative cells, and the general trend 
of telomere shortening, we showed that there is a phase 
transition between a normal state where a cell culture 
becomes senescent to a new state where a cell culture 
can become immortal. 

From our theory, we were able to answer the basic 
questions posed in the introduction. Speciflcally: 

i. We determined the condition for whether a 
cell population ultimately becomes senescent or 
whether it continues to divide ad infinitum. The 
transition between these two regimes is given by 
the condition ej\/ = 1, where em — ADk/Mv^ is a 
dimensionless measure of the relative effect of cell 
division, random fluctuations, and basal loss in the 
length evolution of a telomere. 

ii. We also found that for the case of M = 1 telom- 
ere per cell, the mean time for a cell population to 
become senescent is 



it) 



Xq 



for ei < 1. Here xq is the initial length of the 
telomere, v is the amount by which the telomere 
shrinks by basal loss in each division, and em = 
ADk/Mv^ is a dimension measure of the relative 
importance of cell division to basal loss. 

iii. Finally, we found that the total number of cells pro- 
duced before the entire cell culture becomes senes- 
cent becomes extremely large as e approaches its 
critical value from below, as presented in Eq. (fT4|) . 



Our results may be helpful for understanding how ALT 
cells maintain their telomeres. In these cells, increased 



T-SCE rate, and wide telomere size distribution, and in- 
creased cell lifetimes have been observed jSj] . These obser- 
vations are all natural outcomes from our model. How 
these results depend on the number of short telomeres 
and on the details of the T-SCE process are very impor- 
tant questions that are under investigation. 
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APPENDIX: FIRST-PASSAGE FOR BIASED 
DIFFUSION 

We recall some basic results about first-passage for bi- 
ased diffusion on the positive half line a; > [§] that will 
be used to describe cell proliferation statistics. Accord- 
ing to Eq. ([!]) , the length of each telomere undergoes bi- 
ased diffusion, in the continuum limit, with positive bias 
{v > 0) defined to be directed towards smaller telomere 
length X. An absorbing boundary condition at the origin 
imposes the constraint that when x reaches zero the cell 
effectively becomes senescent and is thus removed from 
the population of dividing cells. 

Let c{x, t) be the concentration of diffusing particles at 
X at time t. The concentration evolves by the convection 
diffusion equation 



dc{x,t) dc{x,t) d^c{x,t) 



dt 



dx 



D- 



dx"^ 



(A.l) 



For the initial condition c(a;,< = 0) = 6{x ~ xo), corre- 
sponding to a single particle starting at xq, the concen- 
tration at any later time is Q 



c{x, t) 



1 



y/AnDt 



^-{x-xo+vty /iDt 



^+vxq/D ^-{x+xo+vty/iDt 



(A.2) 



The second term represents the "image" contribution; 
notice that the bias velocity of the image is in the same 
direction as that of the initial particle. The exponential 
prefactor in the image term ensures that the absorbing 
boundary condition c{x = 0, i) = is always fulfilled. 

From this concentration profile, the first-passage prob- 
ability, namely, the probability for a diffusing particle to 
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hit the origin for the first time at time t, is 



^-{xo-vtY/iDt 



The convective contribution to the flux, — wc, gives no 
contribution because c = at a; = 0. Until the parti- 
cle hits the boundary it stays in the system; hence its 
survival probability is simply 



S{t) = / c(a;, t) dx 
Jo 



erfc 



/ vt -xo \ 



e-o/^erfcf^^) 



(A.4) 



where erfc(z) is the complementary error function. Since 
there is only one absorbing point in the system, the sur- 



vival probability and the first passage probability are re- 
lated by 



Fit) = 



dS(t) 
dt 



(A.5) 



From the asymptotics of the error function [19] , erfc(2:) 
/^/tt z, the long-time behavior of S{t) is given by 



S{t) 



2x0 



'-^ e^^"/2^ g-(.V4D)t^ (A.6) 



As expected intuitively, the survival probability asymp- 
totically decays exponentially with time because the bias 
drives the particle towards the absorbing point. This re- 
sult is used in Eqs. ([6]) & ([7]) to determine the number of 
active cells as a function of time. 
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